Ion irradiation tolerance of graphene as studied by atomistic simulations 
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As impermeable to gas molecules and at the same time transparent to high-energy ions, graphene 
has been suggested as a window material for separating a high- vacuum ion beam system from targets 
kept at ambient conditions. However, accumulation of irradiation-induced damage in the graphene 
membrane may give rise to its mechanical failure. Using atomistic simulations, we demonstrate that 
irradiated graphene even with a high vacancy concentration does not show signs of such instability, 
indicating a considerable robustness of graphene windows. We further show that upper and lower 
estimates for the irradiation damage in graphene can be set using a simple model. 



Since the isolation of graphene in 2004 [l| , a multitude 
of applications have been proposed for this one-atom- 
thick carbon membrane Although graphene is prob- 
ably best known for its unique electronic properties [3|, 
it is also the strongest material ever measured [4] . More- 
over, on the one hand, even a single layer of graphene can 
withstand pressure imposed by a macroscopic amount of 
gas [5] , and does not allow even the smallest atmospheric 
molecules to permeate [6]. On the other hand, porous 
graphene has been considered to be the ultimate mem- 
brane for gas separation f?\ and an ideal material for 
super capasitors [8]. At the same time, graphene is virtu- 
ally transparent to high energy ions, which pass through 
the material without creating substantial damage due 
to negligible interaction cross section [9]. Due to the 
unique combination of strength and impermeability to 
gases, graphene can be used as a window material in ex- 
ternal ion beam experiments with samples which cannot 
be put into vacuum required for the operation of the ion- 
beam system [9], replacing silicon-nitride membranes [10] 
currently used for this purpose. Very recently, a simi- 
lar technology has been demonstrated with graphene ox- 
ide windows for in situ environmental cell photoelectron 
spectroscopy [11]. 

Nevertheless, irradiation damage can accumulate with 
increasing ion dose, so that the operation of graphene 
as an ion-transparent but gas-separating membrane de- 
pends crucially on its ability to withstand continuous ir- 
radiation during the experiment. Although the ion irra- 
diation response of graphene has been recently studied 
both experimentally [12-26] and theoretically [9, 27, 28], 
the atomic-scale details of damage accumulation during 
continuous exposure to the ion beam remain unknown. 
In the experiments, ion doses were either very low creat- 
ing only spatially well-separated defects or so high that 
graphene was completely destroyed, and it was not pos- 
sible to get any insight into damage accumulation pro- 
cess. Similarly, the theoretical work has hitherto concen- 
trated on the effects of individual ion impacts on pris- 
tine graphene. Continuous high-dose irradiation has been 
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taken into account only in a stochastic manner disregard- 
ing the actual dynamics of consecutive ion impacts on a 
defective graphene structure [28]. 

In this Letter, we utilize atomistic simulations to study 
the effects of ion irradiation on graphene with defects 
to understand the details of damage accumulation. We 
demonstrate that irradiated graphene with vacancy con- 
centration of at least 35% does not show any signs of 
structural failure, pointing to considerable stability of 
graphene windows in ion beam experiments. Our results 
can be directly utilized both in estimating the wear of 
graphene windows used to separate ion beam systems 
from volatile targets as well as in designing optimum pa- 
rameters for carving nanopores into graphene membranes 
using a focused ion beam. 

We created our graphene target structures by ran- 
domly removing atoms from the pristine lattice, then 
running a 3 ps annealing simulation at 800 K to allow sat- 
uration of open bonds where possible, and finally relaxing 
the structure to its local energy minimum by quenching 
to K. At this point we also removed all isolated frag- 
ments leaving only the largest continuous atomic network 
to be used as a target for our ion irradiation simulations, 
as illustrated in Fig. [TJ The initial structure consisted of 
1250 carbon atoms. We performed in total 200,000 ir- 
radiation simulations for individually created structures 
with impact points randomly selected in the middle of 
the simulation cell. For each structure, ion species (He, 
Ar, Xe), ion energy {K G [30 eV, 1 MeV]) and initial va- 
cancy concentration (cy) combination, we performed on 
average more than 1050 independent simulations in order 
to collect representative statistics. 

Irradiation simulations were carried out using the 
molecular dynamics (MD) method as implemented in 
the PARCAS simulation code [29] . The carbon-carbon in- 
teraction was described by the reactive bond-order po- 
tential developed by Brenner et al [30|, HH, disregard- 
ing the bond conjugation term which is not important 
for ion irradiation effects [32|. The interactions between 
noble gas ions [33^ and carbon atoms were modeled by 
the universal repulsive potential by Ziegler, Biersack and 
Littmark |3^ . A similar repulsive potential was fitted to 
the carbon-carbon interaction at short distances in or- 
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Figure 1. Overview of the simulations setup for a very high 
vacancy concentration. Bonds between the atoms are drawn 
with sohd hues, (a) An example structure as created with 
the vacancies marked with open squares. Thicker (red) lines 
show isolated islands which were also removed from the tar- 
get structure, (b) The structure after annealing and struc- 
tural optimization. The wide semi-transparent lines at the 
boundaries indicate the heat dissipation area. 



der to properly describe high energy collisions between 
carbon atoms. A few atomic rows at the edges of the 
periodic directions were coupled to the Berendsen ther- 
mostat [355 kept at K in order to model the dissipa- 
tion of heat from the irradiated area^ as indicated in 
Fig. [TJd. We have previously used a similar simulation 
setup for modeling ion irradiation of carbon nanotubes 
Hsl, [33, pristine graphene [9, 27] and hexagonal boron 
nitride mono-layers [38|. After each ion impact we again 
quenched the system into a local energy minimum for 
analysis. 

In Fig. [2] we present sputtering yield (V) for different 
ions and irradiation energies as a function of cy. Note 
that in the limit cy ^ 0, the results approach those for 
pristine graphene [9], as expected. Despite the fact that 
the presented results are averaged over a large number 
of simulations, significant statistical variations remain 
present in the data. However, the qualitative behavior 
of y as a function of vacancy concentration is neverthe- 
less apparent. For instance, it is clear that Y tends to 
decrease with increasing cy due to increased probability 
for the ion to pass through an existing vacancy with an 
energy-dependent slope until the highest considered va- 
cancy concentrations (cy < 35%). Perhaps the most sur- 
prising result is the continuity of all of the curves (except 
for some statistical fluctuations), since one could assume 
that at some point the defective membrane would become 
structurally unstable so that it would break showing an 
abrupt increase in the sputtering yield (corresponding to 
a lost membrane). Our results show that such an insta- 
bility point is not reached within any reasonable vacancy 
concentration. Indeed, a careful analysis of the distri- 
bution of the number of sputtered atoms as a function 
of cy shows no abrupt changes in the relative probabil- 
ities for sputtering smaller or larger numbers of atoms 
at once. The structural stability of the membranes was 
further checked by 1 ns anneal simulations at 1500 K for 
150 structures with different cy after the analysis. No 



apparent instabilities, e.g., disintegration or crumpling, 
were observed. We also point out that although mechan- 
ical properties of perforated graphene are inferior as com- 
pared to the pristine material, the fracture stress remains 
as high as ~ 50 GPa corresponding to nearly 10% strain 
for graphene with up to 2 nm holes [39]. 

An additional surprise is the clearly evident ion energy- 
dependent slope of the data presented in Fig. [2l Based 
purely on geometrical arguments, Y should decrease with 
Cy due to a drop in target density and thus smaller col- 
lision cross section as 

Y{cy) = {l-Cy)Yo, (1) 

where Yq = Y{cy = 0) corresponds to pristine graphene. 
Since this equation does not depend on the energy of 
the impinging ion, it is clear that it cannot completely 
describe the data shown in Fig. [2j Therefore, to quan- 
tify the differences between the data and the geometric 
model, we introduced a new dimensionless variable 7 in 
Eq. [T]to obtain 

Y{cy) = {l-^Cy)Yo. (2) 

Fitting the data presented in Fig. [2] to Eq. [2] allows us 
to systematically analyze the deviation from the geomet- 
ric model as a function of ion energy. The results of the 
fits are presented in Fig. [3l Naturally, the fitted sputter- 
ing yield at zero vacancy concentration (Yq) for each en- 
ergy/ion combination (Fig. [3^) is nearly identical to the 
simulation data for pristine graphene [9|. The modified 
equation reduces to the simple geometric model (Eq. [1]) 
for 7 = 1, whereas negative values indicate increasing 
Y for increasing cy (a positive slope), and 7 = would 
suggest sputtering yield which does not depend on cy. 

As evident from Fig. [3fb), 7 is negative at low ion en- 
ergies. In the case of He and Ar, this occurs only up to 
K = 30 eV, whereas for the heaviest Xe also the 100 eV 
data has a positive slope. For energies immediately above 
these, 7^1 indicating a perfect agreement with the sim- 
ple geometric model. Then, 7 decreases until an apparent 
saturation towards a constant value (within our statisti- 
cal accuracy) at an ion-dependent energy. At the satu- 
ration the average values of 7 are 1.08 ± 0.62, 0.32 ± 0.10 
and 0.53 ± 0.05 for He, Ar and Xe, respectively (notice 
that the statistical uncertainty decreases for increasing 
ion mass). 

The varying values of 7 can be understood taking into 
account how the collision process changes with increas- 
ing ion momentum (and thus energy) [9]. At low en- 
ergies, the collision cross section is large, and the colli- 
sion process is relatively slow. The pre-existing vacan- 
cies in the structure on average lower the binding energy 
of the target atoms which makes an important contri- 
bution to the sputtering when the transferred energies 
are very close to the displacement threshold (22 eV in 
pristine graphene [io', Tlj] and somewhat lower for under- 
coordinated atoms [42]). As our data shows, this effect 
is large enough to overcome the effect of the lowered tar- 
get density at the lowest ion energies. While ion energy 
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Figure 2. Sputtering yield (Y) as a function of vacancy concentration (cv) for (a) He, (b) Ar and (c) Xe ions. Each data point 
is an average Y for all simulations for the same ion species, ion energy and cv- The error bars corresponding to Y show the 
standard deviation and the error bars corresponding to cv are caused by the removal of unconnected fragments after randomly 
removing atoms from pristine graphene. 
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Figure 3. (a) Fitted sputtering yield Yb at zero vacancy con- 
centration (cv = 0) and (b) fitted dimensionless constant 7 
indicating deviation from the simple geometric model. Panel 
(c) presents a magnification of the area marked with a rect- 
angle in panel (b). 



increases, a larger fraction of the recoil atoms will re- 
ceive - during the collision - energies clearly higher than 
the threshold. Therefore, the weaker binding is not sig- 
nificant anymore, and the data agrees better with the 
simple geometric model. However, at even higher ener- 
gies the collision cross section is very small and the time 
during which the ion interacts with the target atoms is 
very short. This means that the target atom will remain 
essentially immobile during the interaction, which corre- 
spondingly becomes symmetric over the graphene plane. 
Therefore, the transferred momentum will almost com- 
pletely be in the in-plane direction. In this case the role 



of the target density is negligible after the initial impact, 
since the displaced target atom will travel in the in-plane 
direction as long as it takes to collide with another target 
atom. While the overall binding of the defective mem- 
brane is lowered, these secondary collisions can displace 
more atoms from the structure than what would hap- 
pen in the case of a pristine target. This sets 7 G [0, 1] 
also for the highest irradiation energies, as can be seen 
in Fig. [3b. Thus, for energies K > 1 keV, the modified 
geometric model can be used to set the upper and lower 
bounds for irradiation-induced damage in graphene. 

We stress that defective graphene layers with rela- 
tively large defects are still hardly permeable [43] for 
small atoms and molecules. Moreover, vacancies in 
graphene tend to partially "heal" themselves by form- 
ing non-hexagonal rings due to bond rotations [Sll, E^ - 
46]. As ion beam can give rise to dissociation of at- 
mospheric molecules, passivation of dangling bonds by 
hydrogen atoms or incorporation of foreign atoms (e.g., 
nitrogen [47]) as substitutional impurities should also de- 
crease the permeability of graphene membranes. 

In conclusion, we have shown that the response of 
graphene to ion irradiation remains consistent to remark- 
ably high vacancy concentrations (up to 35%). Although 
the damaging process varies with the mass of the imping- 
ing ion as well as its energy, we never observed sudden 
breakage of the membrane indicating severe structural in- 
stabilities with respect to the ion irradiation. Naturally, 
whether the mechanical stability of a perforated graphene 
membrane is high enough for a particular application de- 
pends on the actual experimental conditions. However, 
taking this into account, our results can be used to pre- 
dict the effects of ion irradiation on the membrane dur- 
ing the experiment using the presented simple geometric 
model in order to estimate the usability of graphene win- 
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dows. Our results can also be used to design optimal 
conditions for carving nanopores into graphene using a 
focused ion beam. 
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